Shubnikov-de Haas oscillations in YBa 2 Cu 4 08 
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We report the observation of Shubnikov-de Haas oscillations in the underdoped cuprate super- 
conductor YBa2Cu40s (Y124). For field aligned along the c-axis, the frequency of the oscillations 
is 660 ± 30 T, which corresponds to ~ 2.4% of the total area of the first Brillouin zone. The effective 
mass of the quasiparticles on this orbit is measured to be 2.7 ± 0.3 times the free electron mass. 
Both the frequency and mass are comparable to those recently observed for ortho-II YBa2Cu30e.5 
(Y123-II). We show that although small Formi surface pockets may be expected from band structure 
calculations in Y123-II, no such pockets are predicted for Y124. Our results therefore imply that 
these small pockets are a generic feature of the copper oxide plane in underdoped cuprates. 



Understanding the doping evolution of the electronic 
ground state within the superconducting CuC>2 planeš is 
central to the high-T c problem, yet many outstanding is- 
sues remain. Heavily overdoped cuprates, for example, 
behave in manv ways like conventional, highly anisotropic 
metals [H, 0, [Šj], yet quantum oscillations, a character- 
istic signature of a Fermi- liquid ground state jij, have 
never been observed. Moreover, it is usually assumed 
that as hole doping is decreased and the materials evolve 
towards the antiferromagnetic insulating state, the prop- 
erties depart substantially from those of a conventional 
Fermi-liquid. The recent observation [j| of the Shubnikov 
de Haas effect (SdH) in underdoped YBa 2 Cu30 6 .5, with 
T c /T c max ~ 0.6, was therefore a surprising and highly 
significant discovery. 

The frequency of the SdH oscillations in oxygen or- 
dered ortho-II YBa2Cu3C>6.5 (Y123-II) corresponds to a 
Fermi surface pocket with a eross-seetional area (for H||c) 
that is only a small fraetion (~2%) of the total Brillouin 
zone area. The angle dependence of the frequency of 
this oscillation suggests that this pocket is quasi-two di- 
mensional [f|. Density funetional theory (DFT) band- 
strueture calculations [7|, [8J suggest that these pockets 
may be formed from a flat band, originating from the 
CuO chain and BaO layers, which lies very close to the 
Fermi level Ef- However, there are a number of al- 
ternative possibilities such as the formation of pockets 
near the nodal points of the DFT Cu0 2 plane Fermi sur- 
face, where the d-wave superconducting gap is minimal, 
and where the most intense spectral weight is observed 
by angle resolved photoemission spectroscopy (ARPES) 
[9J. These pockets or 'Fermi ares' may occur due to 
promixity to the Mott insulating state [l(| or because 
of reconstruetion of the original Fermi surface resulting 
from the formation of an ordered phase with broken sym- 
metry [ll|, [T3, EH]. Clearly, it is imperative to deter- 
mine whether the SdH signals are unique to Y123-II or 



are an intrinsic feature of the electronic structure of the 
cuprates. 

In this Letter we report the observation of SdH os- 
cillations in another underdoped cuprate superconduetor 
YBa 2 Cu 4 8 (Y124), which has a higher T c (~ 80 K) 
and a higher hole density per planar Cu (p ~ 0.14) than 
Y123-II (p ~ 0.1), in pulsed high magnetic fields up to 
61 Tesla. The frequency of the oscillations is compara- 
ble with those in Y123-II. This similarity, coupled with 
the marked differences in the calculated band-strueture 
of the two materials, strongb/ points to an origin of the 
oscillations which is generic to the underdoped cuprates. 

In contrast to the Y123 family, which has a single CuO 
chain with variable oxygen content, Y124 contains alter- 
nating stacks of Cu02 bilayers and double CuO chains 
that are stable and fully loaded. Single crystalline sam- 
ples (typical dimensions 400 x 80 x 30 /im 3 ) were fiux 
grown in Yg03 erucibles in a partial oxygen pressure 
of 400 bar [14j. Magnetoresistance (MR) measurements 
were performed in the LNCMP pulsed field facility in 
Toulouse. Although several samples were measured, clear 
SdH signals wcre seen in only two samples with the high- 
est T c values (> 80K) and the largest MR response. 
Four-point resistance measurements were taken at a fre- 
quency of 40-80kHz with the current (10-20mA) applied 
I||a (samplc #1, T c ~ 81K) and I||6 (sample #2, T c ~ 
82K). The contact geometry for each crystal is shown 
in the insets of Fig. [1] The voltage aeross the samples 
was recorded using an analogue to digital converter al- 
lowing the data to be post-processed. The samples were 
cooled with either a dilution fridge (#1) or a pumped 
4 He cryostat (#2) housed in two different magnet coils. 
In both cases, the samples were in direct contact with 
the cryogenic liquid at ali temperatures. The magnet 
coils had peak fields of 58T and 61T respectiveb/ and in 
the field range of interest here (45-61T), the maximum 
field rise/fall rates were ~ 2000 T/s. Data for up and 
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FIG. 1: (Color online) Longitudinal and Hali data (p X x, p yy 
and p xy ) for two Y124 samples. The insets show the contact 
configuration for each sample. 



down sweeps showed negligible hysteresis. 

The contact configuration of sample #2 (I||6) allowed 
us to measure both the longitudinal (p yy ) and Hali (p^ y ) 
resistivities by reversing the direction of the applied field, 
whereas for sample #1 (I||ffl), p X x was measured directly. 
In order to improve the signal to noise ratio, data from 
several (typically 2-3) pulses were averaged. Resistivity 
versus field data for both samples are shown in Fig. [TJ 
The resistivities rise rapidly above the irreversibility field 
H m (which is slightly different in the two samples) then 
increase more slowly. Above H lrTl the resistivity is highly 
anisotropic within the plane due to the presence of the 
highly conducting double-chain unit, as found in zero- 
field above T c |l5l | As in Yl 23-11 p^ y is negative and 
varies linearly with H, cxtrapolating to the origin. This 
linear field dependence suggests that the negative p xy is 
a normal state effect. perhaps resulting from a Fermi sur- 
face reconstruction [16J as discussed below, rather than 
a mixed state effect. For both samples and in both the 
longitudinal and transverse components, small period (in 
inverse field) SdH oscillations are visible which are a few 
percent of the total signal at maximum field. 

In Fig.[5]we show the oscillatory component of the MR 
signal for both crystals at various temperatures. Due to 
small differences in noise level for the two field directions, 
our study of the T-dependence of the oscillations in sam- 
ple #2 was carried out with the field in one direction only. 
The resistivity p shown in Fig. [5] for #2 is therefore a 
mbcture of p yy and p*£ . 
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FIG. 2: (Color online). The oscillatory component of the 
magnetoresistance [Ap/ p(H)] for both samples at various 
temperatures. For #1, p = p xx (H) whilst for #2, p is a mix- 
ture of p yy (H) and p^ y (H). In each čase, a linear background 
has been subtracted and the data offset for claritv. The solid 
lines are fits to Eq. {TJ giving m* = 2.7 ± 0.3m e . 



We analyzed the data using the standard theory for 
SdH oscillations in a three-dimensional metal [T3| , which 
is essentially the same as the Lifshitz-Kosevich model for 
the de Haas-van Alphcn oscillations in the magnetisation 
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where A is a constant. The temperature damping factor 
R T = X/{smhX) with X = (2ir 2 k B /he)(m*T/ B) and 
m* is the quasi-particle effective mass. The Dingle fac- 
tor R D = exp(— it/uj c t) — exp[— (nhkF)/{eB£)] deseribes 
the increase in damping as the orbitally averaged mean- 
free-path t decreases or as the average Fermi wavevector 
(&f) of the orbit grows. 

A fit of this equation to the data in Fig.[^]gives a single 
SdH frequency F of 660 ± 30 T for both samples. Us- 
ing the Onsager relation, F = (h/2ire)A, we find that A, 
the eross-seetional area of this orbit, corresponds to 2.4% 
of the average eross-seetional area of the first Brillouin 
zone (27.9kT for H||c). Assuming that the Fermi surface 
or pocket corresponding to this orbit is two-dimensional, 
and that there are n sueh pockets per CuC>2 sheet, Lut- 
tinger's theorem givesp = nabA/An 2 — 0.024(2)n (where 
a, b are the lattice constants). More importantly, A is 
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FIG. 3: (Color online). SdH signal amplitude versus temper- 
ature for sample #2. The solid line is a nt to the standard 
expression for the temperature damping factor, Rt, yielding 
an effective mass of 2.5 ± 0.3m e . The faint dashed lines give 
the approximate range of masses consistent with the present 
data. Inset: Fast Fourier transform of p(l/B,T). 



found to be 1.25±0.07 larger than that observed in Y123- 
II for which p ~ 0.10. Although it is difficult to pinpoint 
precisely the actual number of doped holes p per Cu02 
plane in Y124, taking the average of p fr om three inde- 
pendent estimates (by comparing T c [18| | and the pseu- 
dogap temperature T* UŠI . Hil . [20j ] with Y123 and the 
Seebeck coefficient [2lL l22jV one obtains p ~ 0.14. Thus 
A is found to scale approximately linearly with the nom- 
inal hole densities of Y123-II and Y124. 

As a first approach for determining m*, we neglect the 
field dependence of Rt and fix uj c t/B, F and cf> to the 
values found at the lowest temperature so that changes 
in the values of A represent Rt at an effective field B e = 
55 T. A fit of R T to A(T) is shown in Fig. [3] and gives 
m* = 2.5 ± 0.3m e , where m e is the free electron mass. 
This approach slightly underestimates the mass, and a 
better approach is to fit ali the data at ali temperatures 
to Eq. dXJ) simultaneouslv. These fits, shown by solid lines 
in Fig. [21 give m* — 2.7 ± 0.3m e . This is a heavier mass 
than that observed in Y123-II by a factor 1.4 ± 0.2. 

Fast Fourier transforms of p(l/B,T) are shown in the 
inset to Fig. [31 Both the frequency of the main peak 
and the mass derived from the T-dependent amplitude 
are consistent with the direct fitting deseribed above, al- 
though the former technique is more accurate when so 
few oscillations are present. 

From our fits in Fig. [21 we obtain a Dingle factor 
(e- Q / s , a = 340 ± 100 T" 1 ) for sample #2 which for 
cylindrical pockets gives an average SdH mean-free-path 
£ S dH^ 90 ± 30A. Unfortunateh/, we cannot obtain an in- 
dependent estimate of £t r , the transport mean-free-path 
of the in-plane carriers, as we do not yet know how many 
pockets are present. For sample #1, the magnitude of the 
zero-field residual in-plane resistivity p a o ~ 6 ± 2 pJ!cm, 



estimated by fitting the high-field data to the expression 
p a {H) = p a0 + AH 2 /(1+BH 2 ) that deseribes the form 
of the MR beyond the weak-field limit [SŠj]. Assuming 
there are n pockets per CuC>2 plane (and 2n in total), 
this would give i tI = hc/2nk F e 2 p a0 = 2000/n A (±30%). 
Here, c = 13. 6A is the c-axis lattice spacing and kp = 
1.4 nm -1 , as determined from the SdH frequency. 

The above analysis implies that there is either a large 
number of Fermi pockets in Y124, or that the Dingle 
factor grossly underestimates the transport lifctime [2J| . 
The latter could be attributed either to the stronger effect 
of small-angle scattering on the dephasing of quantum os- 
cillations, as seen in simple metals 17j, or to additional 
damping in the supercondueting mixed state. Quantum 
oscillations have been observed in the supercondueting 
state of a number of different materials (see for exam- 
ple, Refs. 2Š], 26|). Whilst the frequencies and effective 



masses are found to be the same as in the normal state, 
an additional damping is observed which is related to the 
magnitude of the field dependent energy gap. We note 
however, that if the orbits observed here originated from 
the nodal regions of the Fermi surface, as suggested by 
ARPES [9J] , then it may be expected that this additional 
damping is small because of the d-wave symmetry of the 
supercondueting gap funetion. 

As mentioned above, conventional DFT band- 
strueture calculations 0, HI for Yl 23-11 show that with 
small rigid shifts of bands (AE ~ ±20 meV), a flat band 
arising from the CuO chain layer and apical oxygen (BaO 
layer) can break through Ep giving rise to small pockets, 
with frequencies rising up to ~ 500 T and bare masses 
~ 1 — 3 m e . To investigate the possibility of similar pock- 
ets arising in Y124 we calculated the band-strueture of 
Y124 using the Wien2K package which is an imple- 
mentation of a full-potential, augmented plane wave plus 
local orbital scheme. The crvstal strueture from Ref. 2Š|, 
and 10 4 fc-points in the full Brillouin zone were used for 
the consistency cycle. 

The band-strueture and Fermi surface are shown in 
Fig. [H and are essentially the same as those reported pre- 
viously [2^, [Š3| ■ The Fermi surface consists of two large 
hole-like quasi-2D cylinders centered on the zone corner 
(J) with mostly plane character and two warped quasi- 
1D sheets with mostly chain character. Cruciallv, the 
position of the above mentioned CuO/BaO band is now 
~ 400 meV below Ep (at the J point), and so a very large 
shift of the energy of this band would be needed to pro- 
duce Fermi surface pockets of similar origin to those cal- 
culated for Y123-II. In fact, Ep needs to be shifted down 
by ~ 80meV before any new closed orbits are formed for 
H || c. The orbit formed is centered around the zone cen- 
ter (r) (see Fig. |4]i) and has a high frequency (~3500 T 
for AEf = —100 meV) so cannot feasibly be responsiblc 
for our observations. Given the similarity of our SdH rc- 
sults for Y124 and those for Y123-II, it therefore seems 
unlikely that the oscillations in either material can be 
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FIG. 4: (Color online). a) Band-structure of Y124 along sym- 
metry directions in the basal plane. The color scale reflects 
the band character: C=CuO chain character (including api- 
cal oxygen) P=Cu02 plane character. b) Brillouin zone of the 
Ammin space group of Y124, showing the symmetry points 
labelled in a), c) Fermi surface of Y124 viewed along the c- 
axis. The shaded region is the first Brillouin zone. d) As for 
c) but with Ef shifted down by 100 meV. 



explained by conventional band-structure theory. 

As proposed for Yl 23-11 5], an alternative to the con- 
ventional band-structure approach is to associate our 
SdH orbits with Fermi surface reconstruction, due, e.g. 
to spin-density-wave [ll|, [H or d-density -wave [12fl for- 
mation. Unfortunately, it is not possible at this stage 
to have a complete picture of the (reconstructed) Fermi 
surface until we know precisely how many SdH orbits are 
present in Y124. In this regard, measurements to highcr 
fields should prove very informative. 

In summary, we have observed Shubnikov-de Haas os- 
cillations in the underdoped cuprate Y124. The fre- 
quency of the oscillations is comparable with those ob- 
served in Yl 23-11, scaling approximately by the ratio of 
hole doping levels in the two compounds. It is likely 
that these oscillations are a generic feature of underdoped 
cuprates perhaps resulting from a reconstruction of the 
CuC>2 plane Fermi surface, although as both materials in 
which SdH oscillations have been observed so far (Y123- 
II and Y124) contain quasi-lD CuO chain structures, an 
alternative origin cannot be completely ruled out at this 
tirne. Nevertheless our results provide precise informa- 
tion on how the pocket area and effective masses vary 
as a function of doping, providing a strong constraint on 
theories which seek to explain their electronic structure. 
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